def proba_binom(n,p,k) :
    coeff_bino = pascal(...)[k]
    proba = coeff_bino*...**(k)*(1-p)**(...)
    return proba